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ABSTRACT 

We investigate the peculiar velocities predicted for galaxy clusters by theories in the 
cold dark matter family. A widely used hypothesis identifies rich clusters with high 
peaks of a suitably smoothed version of the linear density fluctuation field. Their 
peculiar velocities are then obtained by extrapolating the similarly smoothed linear 
peculiar velocities at the positions of these peaks. We test these ideas using large high 
resolution N-body simulations carried out within the Virgo supercomputing consor- 
tium. We find that at early times the barycentre of the material which ends up in a rich 
cluster is generally very close to a high peak of the initial density field. Furthermore the 
mean peculiar velocity of this material agrees well with the linear value at the peak. 
The late-time growth of peculiar velocities is, however, systematically underestimated 
by linear theory. At the time clusters are identified we find their rms peculiar velocity 
to be about 40% larger than predicted. Nonlinear effects are particularly important in 
superclusters. These systematics must be borne in mind when using cluster peculiar 
velocities to estimate the parameter combination crgfi'^'®. 

Key words: cosmology: theory - large-scale structure of Universe, cosmology: theory 
- dark matter, galaxies: clusters 



1 INTRODUCTION 

The motions of galaxy clusters are thought to result from 
gravitational forces acting over the very large scales on which 
superclusters are assembled. The rms deviations from unifor- 
mity on such scales appear to be small, and so may be ade- 
quately described by the linear theory of fluctuation growth. 
For a linear density field of given power spectrum the rms 
peculiar velocity is proportional to agOP'^ where Q is the cos- 
mic density parameter and as , the rms mass fluctuation in a 
sphere of radius 8 /i~^Mpc, is a conventional measure of the 
amplitude of fluctuations (e.g. Peebles 1993). (As usual the 
Hubble constant is expressed as Hq = 100 /i km/sec/Mpc.) 
Distance indicators such as the Tully-Fisher or D„-(7 rela- 
tions allow the peculiar velocities of clusters to be measured, 
thus providing a direct estimate of this parameter combina- 
tion (see, for example, Strauss & Willick 1996). 

Essentially the same parameter combination can also be 
estimated from the abundance of galaxy clusters (e.g. White 
et al. 1993) and a comparison of the two estimates could in 
principle provide a check on the shape of the assumed power 
spectrum and on the assumption that the initial density fleld 



had gaussian statistics. In practice this is difficult because 
of the uncertainties in relating observed cluster samples to 
the objects for which quantities are calculated in linear the- 
ory or measured from N-body simulations. The standard 
linear model was introduced by Bardeen et al. (1986; here- 
after BBKS). It assumes that clusters can be identified with 
"sufficiently" high peaks of the linear density field after con- 
volution with a "suitable" smoothing kernel. The peculiar 
velocity of a cluster is identified with the linear peculiar ve- 
locity of the corresponding peak extrapolated to the present 
day. In the present paper we study the limitations both of 
this model and of direct N-body simulations by comparing 
their predictions for clusters on a case by case basis. 

In the next section we summarize both the linear pre- 
dictions for the growth of peculiar velocities and the BBKS 
formulae for the values expected at peaks of the smoothed 
density field. Section 3 then presents our set of N-body sim- 
ulations and outlines our procedures for identifying peaks 
in the initial conditions and clusters sd, z — 0. Section 4 be- 
gins by studying how well peaks correspond to the initial 
barycentres of clusters; we then show that the smoothed 
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linear velocity at a peak agrees well with the mean linear 
velocity of its cluster; finally we show that the growth of 
cluster peculiar velocities is systematically stronger at late 
times than linear theory predicts. A final section presents a 
brief discussion of these results. 



elliptic integrals. He also shows that the above approxima- 
tions always have fractional errors better than 2% if > 0.1. 
We therefore work with the simpler approximate formulae 
in the present paper. 



2.2 The Velocities of Peaks 



2 LINEAR PREDICTIONS FOR THE 
PECULIAR VELOCITIES OF PEAKS 

2.1 The Growth of Peculiar Velocities 

According to th e linear theor y of gravitational instability in 
a dust universe (Peebles 1993), the peculiar velocity of every 



mass element grows with cosmic scale factor a as 

V oc aD , (1) 

where a{t) is obtained from the Friedman equation 

2 



fio a ^ + (1 - fio - Ao) a ^ + Ao , 



(2) 



and D{t) is the growth factor for linear density perturba- 
tions, 5{x,t) = Z)(t)(5o(x). Qo and Ao are the density pa- 
rameter and the cosmological constant at z = 0, respectively, 
and we define a = 1 at this time. A number of accurate ap- 
proximate forms are known for the relations between D and 
a and can be used to cast the scaling of eq. into a more 
convenient form. We write 
• ^ dL» _ dD da 

dt ' da dt' ^■^> 

and substitute for da/dt from the Friedman equation 
Lahav et al. (1991) give an approximation for dD/da in the 
combination 

0.6 



/(«) 



dD a 
'daD 



fio a-'-^ + (1 - - Ao) a-^ + Ao 



■(4) 



For a = 1 this gives the standard factor / ~ Qq'^' 
appears when predicting the peculiar velocities produced by 
a given overdensity field. Carroll et al. (1992) used this result 
to derive an approximation for D{a) itself. 



D ^ ag{a) , 
where 



2 n^/^ia) ~ A(a) + + + 



with 

A{a) 



a + Qo{l - a) + A(a3 - a) 



a + Oo(l — a) + A(a^ — a) 



(5) 
(6) 

(7) 

(8) 



Combining these equations, we obtain an explicit approxi- 
mation for the growth of peculiar velocities. 



V oc /(a) g{a) a^-^Sloa ^ + (1 — fio — Ao) a^^ 4- Ao . 



(9) 



For the simple Einstein-de Sitter case where flo = 1 and 
A = 0, these formulae reduce to the exact results D = a (x 
1?^^ and V oc ^Ja. 

Recently, Eisenstein (1997) has shown that the exact 
solutions for D and /(a) can be given explicitly in terms of 



The idea that the statistical properties of nonlinear objects 
like galaxy clusters can be inferred from the initial linear 
density field was developed in considerable detail in the mon- 
umental paper of Bardeen et al. (1986; BBKS). If the initial 
fiuctuations are assumed to be a gaussian random field, they 
are specified completely by their power spectrum, P(A;). Sim- 
ilarly, any smoothed version of this initial field is specified 
completely by its own power spectrum, P(fc)Vy^(fci?), where 
W(kR) is the Fourier transform of the spherical smoothing 
kernel and iZ is a measure of its characteristic radius. In par- 
ticular, BBKS showed how the abundance and rms peculiar 
velocity of peaks of given height can be expressed in terms 
of integrals over P(k)W'^(kR). The difiiculty in connecting 
this model with real clusters lies in the ambiguity in deciding 
what smoothing kernel, characteristic scale, and peak height 
are appropriate. Typically the smoothing kernel is taken to 
be a gaussian or a top-hat, R is chosen so that the kernel 
contains a mass similar to the minimum mass of the cluster 
sample, and the height is assumed sufficient for a spherical 
perturbation to collapse by z = 0. 

The smoothed initial peculiar velocity field is isotropic 
and gaussian with a three-dimensional dispersion given by 



(10) 



where, in the notation of BBKS, aj is defined for any integer 
j by 



^'iiR) = TT-T I Pik)W^(kR)k^'+^dk. 



(11) 



The rms peculiar velocity at peaks of the smoothed density 
field difi^ers systematically from (t„; BBKS show that it is 
given by 



a^(R) = a4R)^l-a^/alal,. 



(12) 



Note that this expression does not depend on the height of 
the peaks. As shown in BBKS, the velocities of peaks are 
statistically independent of their height. 

Throughout this paper we will approximate the power 
spectra of CDM models by the parametric expression of 
Bond & Efstathiou (1984), 



P{k,r) 



Ak 



{1 + [ak/r + (6fc/r)^/2 + (cfc/r)2]-}2/'' 



(13) 



where a = (6.4/r) /i"^ Mpc, b = (3.0/r) ft'^ Mpc, c = 
(1.7/r) /i^^ Mpc, u — 1.13, and the shape parameter F is 
given for the models discussed below by 

f noV[0.861+3.8(m?ord)2''=']'''2 for rCDM, 
1 Qoh otherwise. 

In the rCDM case mio is the r-neutrino mass in units of 
10 keV and Td is its lifetime in years (White et al. 1995) . A 
detailed investigation of this model can be found in Bharad- 
waj & Sethi (1998). For the cosmologies used here, we take 
the values shown in Table 1. Detailed calculations of the 
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power spectrum are actually better fit by slightly smaller 
values of F than we assume (Sugiyama 1995). 

The normalisation constant in equation ( |l3| ) can be re- 
lated to the conventional normalisation erg by noting that 
as = cro(8/i~^Mpc) and using equation ( pl| ) with a top-hat 
window function, Wth{x) = 3(xsina: — cos3;)/a;'^. This cor- 
responds to the linear fluctuation amplitude extrapolated to 
z = and can be matched to observation by fitting either to 
the cosmic microwave background fluctuations measured by 
COBE or to the observed abundance of rich galaxy clusters. 
The models of this paper are normalised using the second 
method (c.f. Eke et al. 1996), as reflected by the as values 
given in Table 1 together with the other parameters defining 
the models. 

In the following linear density fields are smoothed ei- 
ther with a top-hat or with a gaussian. In the latter case 
the window function is Wg{x) = exp[— a::^/2]. It is un- 
clear for either filter how R should be chosen in order to 
optimize the correspondance between peaks and clusters. 
We follow previous practice in assuming that cluster sam- 
ples contain all objects with mass exceeding some thresh- 
old Mmin, and then choosing R so that the filter contains 
Mmin- Hence Mmin ~ 4-jTpR^/S in the top-hat case and 
Mmin = {2-Kf''^pR^ in the gaussian case. The simulations 
analysed here have f2o = 0.3 or 1.0, and we will isolate clus- 
ter samples limited at Mmin ~ 3.5 x lO^*/i~^M0, the value 
appropriate for Abell clusters of richness one and greater 
(e.g. White et al. 1993). A detailed discussion of filter- 
ing schemes can be found in Monaco (1998) and references 
therein. 

Table 2 gives characteristic filter radii R and values of 
av and a^ from equations (11) and (13) for both smoothings 
and for all the cosmological models we consider in this pa- 
per; the velocity dispersions are extrapolated to the linear 
values predicted at z = 0. The difference between a^ and a^ 
has often been ignored in the literature when predicting the 
peculiar velocities of galaxy clusters (e.g. Croft & Efstathiou 
1994; Bahcall & Oh 1996; Borgani et al. 1997); for our mod- 
els the two differ by about 15%. Notice also that with our 
choice of filter radii, gaussian smoothing predicts rms pecu- 
liar velocities about 10% smaller than top-hat smoothing. 



3 THE SIMULATIONS 
3.1 The Code 

The Virgo Consortium was formed in order to study the 
evolution of structure and the formation of galaxi es using 



Table 1. The Virgo models 



the latest generation of parallel supercomputers (Jenkini 



et al. 1996). The code used for the simulations of this pa- 
per is called HYDRA. The original serial code was developed 
by Couchman et al. (1995) and was parallelized for CRAY 
T3D's by Pearce et al. (1995) (A detailed description can be 
found in Pearce & Couchman (1997). T3D-HYDRA is a par- 
allel adaptive particle-particle/particle-mesh (AP^M) code 
implemented in CRAFT, a directive-based parallel Fortran 
developed by CRAY. It supplements the standard P^M algo- 
rithm (Efstathiou et al. 1985) by recursively placing higher 
resolution meshes, refinements, over heavily clustered re- 
gions. Refinements containing more than ~ 10^ particles are 
executed in parallel by all processors; smaller refinements 
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are completed using a task farm approach. This T3D ver- 
sion currently includes an SPH treatment of gas dynamics, 
but this was not used for the simulations of this paper. 

A second version of HYDRA, based on CRAY's shared 
memory and message passing architecture, has been written 
by MacFarland et al. (1997). This can run on CRAY T3E's 
but does not currently include refinement placing. 

The simulations used here were run on the Cray T3D 
and T3E supercomputers at the computer center of the Max 
Planck Society in Garching and at the Edinburgh Parallel 
Computing Centre. 



3.2 The Simulation Set 

A set of four matched N-body simulations of CDM universes 
was completed in early 1997. Each follows the evolution of 
structure within a cubic region 240ft~^Mpc on a side us- 
ing 256^ equal mass particles and a gravitational softening 
of 30/i~^kpc. The choices of cosmological parameters corre- 
spond to standard CDM (SCDM), to an Einstein-de Sitter 
model with an additional relativistic component (rCDM), 
to an open CDM model (OCDM), and to a flat low density 
model with a cosmological constant (ACDM). A list of the 
parameters defining these models is given in Table 1. 

In all models the initial fiuctuation amplitude, and so 
the value of as, was set by requiring that the models should 
reproduce the observed abundance of rich clusters. Further 
details of this choice and of other aspects of the simulations 
can be found in Jenkins et al. (1998). Note that each Fourier 
component of the initial fiuctuation field had the same phase 
in each of these four simulations. As a result there is an 
almost perfect correspondance between the clusters in the 
four models. 

Because of their finite volume, these simulations contain 
no power at wavelengths longer than 240/i~^Mpc. Further- 
more, Fourier space is sampled quite coarsely on the largest 
scales for which they do contain power, and so realisation 
to realisation fluctuations on these scales can be significant. 
The size of the effects can be judged from Table 2 where we 
list the values of a^ and ap obtained for each model when 
the theoretical power spectrum is replaced in equations (11) 
and (13) by the initial power spectrum of the model itself. 
These are systematically smaller than the values found be- 
fore. The difference is primarily a reflection of the loss of 
large-scale power. 



3.3 The Selection of Peaks 

We identify peaks in the initial conditions of the simula- 
tions by binning up the initial particle distribution on a 
128^ mesh using a cloud-in-cell (CIC) assignment and then 
smoothing with a gaussian or a top-hat with characteristic 
scale R corresponding to Mmin ~ 3.5 x lO^''/i~^M0. A peak 
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Table 2. For each of the models, the following quantities are given: the radius R (second and fifth column) of the filter used in eq. (12); 
the three-dimensional velocity dispersions <7v and dp (third, fourth, sixth, and seventh column) obtained using eq.s (11) and (13) with 
the given filter radii; the three— dimensional velocity dispersions (Jy and <7p (eighth, ninth, eleventh, and twelvth column) obtained using 
eq.s (11) and (13) with the given filter radii and the power spectra of the simulations themselves; the rms linear overdensity A (tenth 
and thirteenth column) smoothed with the given filter radii and extrapolated to z = 0; the number of clusters Nqi (fourteenth column) 
found in the simulations at 2 = 0; the three-dimensional velocity dispersions of peaks (fifteenth and sixteenth column) in the initial 
conditions of the simulations using the given filters; the three-dimensional linear velocity dispersions of clusters extrapolated to z = 0; 
and the three-dimensional measured velocity dispersion of clusters at 2: = 0. The radii are given in Mpc//i, the velocity dispersions in 
km/sec. Top Hat and Gaussian filters are abbreviated as TH and G, respectively. 
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Figure 1. The mass of the clusters in our simulations against the height of the corresponding peaks in the initial conditions, once these 
are smoothed with a top-hat with the characteristic radius listed in Table 2. All clusters with mass greater than 3.5 x 10^* Mq and 
all peaks with height greater than u = 1.5 are shown. There are 351, 239, 84, and 83 unmatched peaJcs in the SCDM, rCDM, ACDM, 
and OCDM model, respectively. 



is then taken to be any grid point at which the smoothed 

density is greater than that of its 26 nearest neighbours. 
The dimensionless height of a peak, i^, is defined by divid- 
ing its overdensity by the rms overdensity. A, which we list 
in Table 2. Again, within the matched set there is a close 
correspondance between the peaks found in the four mod- 
els. In addition, the peaks found with gaussian smoothing 
correspond closely to those found with top-hat smoothing. 

Particle pecuUar velocities are binned up and smoothed 
in an identical way and the peculiar velocity of a peak is 



taken to be the value at the corresponding grid point. In 
Table 2 we list the rms peculiar velocity of the peaks found 
in each model. Again this is scaled up to the value expected 
at 2 = according to linear theory. It differs slightly from 
the value predicted by inserting the power spectrum of the 
simulation directly into equation (13) because there are re- 
alisation to realisation fluctuations depending on the phases 
of the Fourier components. As it should, the rms peculiar 
velocity averaged over all grid points agrees very well with 
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the value found by putting the simulation power spectrum 
into equation (11). 

3.4 The Selection of Clusters 

We define clusters in our simulations in the same way as 
White et al. (1993). High-density regions at z = are located 
using a fricnds-of-frionds group finder with a small linking 
length (b=0.05), and their baxycentres are considered as 
candidate cluster centres. Any candidate centre for which 
the mass within 1.5 Mpc exceeds Mmin is identified as a 
candidate cluster. The final cluster list is obtained by delet- 
ing the lower mass candidate in all pairs separated by less 
than 1.5 Mpc. In the following we will normally consider 
only clusters more massive than Mmin = 3.5 x 10^*/i~^Mo. 
The number of clusters found in each simulation is listed 
in Table 2. As already noted, the individual clusters in the 
different simulations of the matched set correspond closely. 
Despite the normalisation to cluster abundance it appears 
as though the SCDM model has significantly more clusters 
than the others. This is a reflection of its steeper power 
spectrum together with the value of Mmin we have chosen. 
For Mmin = 5.5 x lO^*/i~^M0 all the models have about 20 
clusters. 

We define the peculiar velocity of each cluster at 2 = 
to be the mean peculiar velocity of all the particles within 
the 1.5 Mpc sphere. The peculiar velocity of the cluster 
at earlier times is taken to be the mean peculiar velocity of 
these particles. Consistent with this, we define the position 
of the cluster at each time to be the barycentre of this set of 
particles. At 2 = this is very close to, but not identical with 
the cluster centre as defined above. We give the rms values 
of the initial (linear) and final {z = 0) peculiar velocities 
of the clusters in each of our models in Table 2. The initial 
values have been scaled up to the linear values predicted 
at = 0. It is clear that these substantially underestimate 
the actual values, a result we discuss in more detail below. 
We note that the present-day properties of clusters in these 
simulations are considered in much more detail in Thomas 
et al. (1998). 



4 COMPARISON OF THE PEAK MODEL 
WITH SIMULATIONS 

4.1 The Cluster-Peak Connection 

The extent to which dark haloes can be associated with 

peaks of the smoothed initial density field is somewhat con- 
troversial. Frenk et al. (1988) concluded that, for appropri- 
ate choices of filter scale and peak height, the correspon- 
dance is good, whereas Katz et al. (1993) claimed that 
"there are many groups of high mass that are not associ- 
ated with any peak" . The result of correlating the peaks in 
the initial conditions of our simulations with the initial po- 
sitions of our clusters is illustrated in Fig. 1. We consider 
a peak and a cluster to be associated if their separation is 
less than 4/i~^ Mpc (comoving). We find that the barycen- 
tcrs of 70% and 80% of the clusters with masses exceeding 
3.5 X lO'^'* Mq lie within of a peak with v > 1.5 for the 
low and high models, respectively. 

Fig. 1 shows that there is, as expected, a correlation 



between the height of a peak and the mass of the corre- 
sponding cluster. In addition, combining the peak heights 
with the A values from Table 2, we see that the extrap- 
olated linear overdensities of the peaks at redshift zero are 
similar but somewhat larger than the threshold value of 1.69 
used in the standard Press-Schechter approach to analysing 
structure formation. 



4.2 Linear Peculiar Velocities of Peaks and 
Clusters 

Given the excellent correspondance between peaks of the 
smoothed linear density field and the initial positions of clus- 
ters, it is natural to compare the smoothed peculiar velocity 
at a peak with the mean initial peculiar velocity of its asso- 
ciated cluster. We show such a comparison in Fig. 2, again 
based on top-hat smoothing of both position and peculiar 
velocity fields using the characteristic radii listed in Table 2. 
All velocities arc scaled up to the expected value at 2 = ac- 
cording to linear theory. The correlation is clearly excellent 
in all cases, and is similar if gaussian rather than top-hat 
smoothing is used. The rms difference in peculiar velocity 
between a cluster and its associated peak is 16%, 16%, 23%, 
and 17% of the corresponding ap value listed in Table 2 
for the OCDM, ACDM, SCDM and tCDM simulations re- 
spectively. The somewhat larger percentage for the SCDM 
model is probably a consequence of the greater influence of 
small-scale power in this case. 

4.3 The Growth of Cluster Peculiar Velocities 

If cluster peculiar velocities grew according to linear theory 
the scaled initial velocities discussed in the last section and 
plotted in Fig. 2 would correspond to the actual velocities 
of the clusters at 2 = 0. In Fig. 3 we show scatter diagrams 
in which these two velocities are plotted against each other. 
It is evident that in fact the agreement is quite poor and 
that there is a systematic trend for the true cluster velocity 
to be larger than the extrapolated linear value. This is re- 
flected in the substantial difference between the rrns values 
of these two quantities listed in Table 2. It is presumably 
a consequence of nonlinear gravitational forces accelerating 
the clusters. 

Some confirmation of this is provided by Fig. 4 where we 
plot the peculiar velocity in units of its initial value for five 
clusters from each of our cosmologies. At early times the 
peculiar velocities all grow as expected from linear theory 
(indicated in the figures by a dotted line) but at later times 
the behaviour is more erratic and most clusters finish with 
larger velocities than predicted. 

Further evidence that late-time nonlineax effects are re- 
sponsible for this discrepancy comes from Fig. 3. In this 
plot all clusters that have a neighbour within 10/i~^Mpc 
arc indicated with a diamond while more isolated clusters 
are indicated by a cross. It is evident that deviations from 
linear theory are substantially larger for the "supercluster" 
objects than for the rest. These objects also have systemat- 
ically larger peculiar velocities at 2 = 0. Their rms peculiar 
velocity is around 20 to 30% larger than that of the sample 
as a whole. 

For the tCDM model, we have run a second realization 
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Figure 2. The initial peculiar velocities of clusters in each of our four cosmogonies are compared to the linear peculiar velocities of their 
associated peaJjs. The linear peculiar velocity field was smoothed with a top-hat in the same way as the density field in order to obtain 
the peak peculiar velocities. 



of the power spectrum. We have extracted a cluster sample 
in the same fashion as described above. The rms pecuUar ve- 
locity of the clusters at z = is 0-2=0 = 511 km/sec. The ex- 
trapolated rms linear peculiar velocity is (Tz=o = 394 km/sec. 
These numbers are very close to the values obtained for the 
first realization. Although two simulations are not a good 
statistical sample, we conclude that there is no realization 
dependence of the mis-match between the extrapolated lin- 
ear and the actual peculiar velocities of galaxy clusters. 

It might be thought that this anomalous acceleration 
of clusters at late times was a consequence of the relatively 
small radius, 1.5/i~^Mpc, which we use to define our clusters. 
Material could, perhaps be ejected asymmetrically from this 
region during the merging events by which clusters form. We 
have searched for such effects by redefining clusters to be all 
the material contained within a radius of 3 or 5 Mpc 
and then repeating the analysis for the same set of objects 
as before. In most cases this turned out to make very little 
difference to either the initial or the final velocities mea- 
sured, and it did nothing to reduce the discrepancy between 
them. The relevant nonlinear effects arc acting on signifi- 
cantly larger scales. Wo repeated this procedure going as far 
out as 25 /t^^ Mpc from the cluster center. At a radius of 
10 Mpc, the difference between the rms peculiar veloc- 
ity and the extrapolated rms linear peculiar velocity is only 
10%. By a radius of 20 Mpc, the numbers have finally 
converged. 

The discrepancy between the rms peculiar velocity of 
clusters and their extrapolated rms linear peculiar velocity is 



independent of any smooting of the density field. With our 
choice of smoothing filter, the linear peculiar velocities of 

our clusters match those of their associated peaks as well as 
the rms value predicted by linear theory when the simulated 
realization of the power spectrum and the proper expression 
for the peculiar velocities (eq. 13) is used. Previous work 
(e.g. Borgani et al. 1997) has tried to match N-body data 
with linear theory by tuning the filter scale. Our results 
undermine the physical basis for such procedure. 



5 CONCLUSIONS 

We have investigated the peculiar velocities predicted for 
galaxy clusters by theories in the Cold Dark Matter family. 
A widely used hypothesis identifies rich clusters with high 
peaks of a smoothed version of the linear density fiuctuation 
field. Their peculiar velocities are then obtained by extrapo- 
lating the similarly smoothed linear peculiar velocities at the 
positions of these peaks. We have tested this using a set of 
four large high-resolution N-body simulations. We identify 
galaxy clusters at z = and then trace the particles they 
consist of back to earlier times. In the initial density field, 
the barycenters of 70% and 80% of the clusters with masses 
exceeding 3.5 x 10^^ Mq lie within 4/i~^ Mpc (comov- 
ing) of a peak with u > 1.5 for the low and high O models, 
respectively. Furthermore, the mean linear peculiar velocity 
of the material which forms a cluster at ^; = agrees well 
with the value at that peak. 
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Figure 3. The initial peculiar velocities of clusters in each of our four cosmogonies, scaled up to z = using linear theory, are compared 
to their actual peculiar velocities at z = 0. Diamonds denote clusters which have a neighbour within 10/i~^Mpc while crosses denote 
more isolated clusters. 



However, the late-time growth of peculiar velocities is 
systematically underestimated by linear theory. At the time 
clusters are identified, i.e. at 2: = 0, we find that the rms pe- 
culiar velocity is about 40% larger than predicted. Nonlinear 
effects are particularly important in superclusters; the rms 
values for clusters which are members of superclusters are 
about 20% to 30% larger than those for isolated clusters. 
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